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The Lindblad equation is an evolution equation for the density matrix in quantum theory. It 
is the general linear, Markovian, form which ensures that the density matrix is hermitian, trace 1, 
positive and completely positive. Some elementary examples of the Lindblad equation are given. The 
derivation of the Lindblad equation presented here is "simple" in that all it uses is the expression 
of a hermitian matrix in terms of its orthonormal eigenvectors and real eigenvalues. Thus, it is 
appropriate for students who have learned the algebra of quantum theory. Where helpful, arguments 
are first given in a two-dimensional hilbert space. 

PACS numbers: 



I. INTRODUCTION 

The density matrix p is a useful operator for quantum 
mechanical calculations. For a given system, one may 
be unsure about what is the state vector. If the pos- 
sible state vectors and their associated probabilities are 
{\'^i)iPi}, one creates the proper^ density matrix 

It is hermitian: = p. It is trace 1: 1 = Trp = 
J2j{'t'j\p\4'j)^ where {I'/'j)} are an arbitrary complete or- 
thonormal set of vectors. It is positive: {v\p\v) > for 
an arbitrary vector \v). All these properties can easily be 
verified from Eq. ([!]). 

One can use the density matrix to conveniently calcu- 
late probabilities or mean values. If a measurement is set 
up to result in one of the eigenstates of an operator 
O, that outcome's probability is and the mean 

eigenvalue of O is TrO/o. 

Because each individual state vector evolves unitarily 
under the system Hamiltonian H (assumed for simplic- 
ity here to be time independent), [ip, t) = exp —ifH]^, 0), 
the density matrix in Eq. ([T]) satisfies the evolution equa- 
tion 

j^p{t) = -i[H, p{t)] ^ Pit) = e-^HV(0)e^«*. (2) 

The operator acting on p{0) is often called a superoper- 
ator^ since it describes a linear transformation on an 
operator: it operates on both sides of p, so to speak. 

The case sometime arises where the system S under 
consideration is a subsystem of a larger system S + S' , 
and 5" is not measured. The pure (so-called because it 
is formed from a single state vector) density matrix for 
the joint system is 

R(i) = E (*) I Ix™) E (*) * ('^^^ I I , 
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where {I'/'i)} {|Xm)} are orthonormal bases for S, S' re- 
spectively, and J2im |C'im(i)P = 1- By taking the trace 
of R with respect to 5", one arrives at an improper^ 
density matrix for S from which predictions can be ex- 
tracted: 

Pit) =Y,C^Ut)CJ,n{tr\^^){<|>J\■ 
ijm 

One can easily see that this is hermitian, trace 1 and 
positive. However, while the density matrix of 5 S" 
evolves unitarily, the density matrix of the subsystem 
S evolving under the influence of S" generally does not 
evolve unitarily. Nonetheless, sometimes dpit) /dt can be 
written in terms of p for a range of times earlier than t. 
Sometimes that range is short compared to the time scale 
of evolution of p so that one may make an approximation 
whereby dpit)/dt depends linearly just on pit). This is 
quite useful, and is what shall be considered in this paper. 

In this case, the evolution equation is highly con- 
strained by the requirements on pit), to be satisfled at 
all times: hermiticity, trace 1 and positivity. (The latter 
proves too general to simply implement, so a stronger 
requirement is imposed, called complete positivity — 
see Section |IV[ ).) The result, for an TV-dimensional 
Hilbert space, is the Lindblad[5] (or Lindblad- Gorini- 
Kossakowsky-Sudarshan|4l) evolution equation for the 
density matrix: 

J^pit)^-^[^l,pit)] 
-- ^ [L"tL"p(i) -f- p(t)L"tL" - 2L"p(i)L"t]. 

Q = l 

(3) 

In Eq.([3|, the hamiltonian H is an arbitrary hermitian 
operator, but the Lindblad operators {L"} are com- 
pletely arbitrary operators. 

Actually, as shall be shown, there need be no limitation 
on the number of terms in the sum in Eq.([3|, but this 
can always be reduced to a sum of N'^ — 1 terms. 

It is not a necessary condition, but if the equation is 
to be time-translation-invariant, the operators are time- 
independent. 
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II. LINDBLAD EXAMPLES 

Before deriving Eq.(|3|, we give a few examples of the 
non- unitary evolutions it describes. Since unitary evolu- 
tion is well known, we shall let H = 0. It shall be seen 
that relaxation to some equilibrium (constant) density 
matrix is readily described. 

For simplicity, four of the five examples shall be in a 
N — 2 hilbert space (a restriction that is readily lifted). 



A. Random Phases 

Consider a state vector written in a basis whose phase 
factors undergo random walk. 

Given an initial state vector jV') = a\(j)i) + b\<j)2) (|ap + 
|6p = 1), suppose at time t, it has evolved to 



with probability 



P{ei,e2)d0ide2 



1 /)2 1 n2 



The density matrix is 

poo poo 

p{t) = / / de^d92P{eue2m{i'\ 



J —oo J — oo 

= ian0i)(0ii + i6n<^2)(</'2i 

+e-5*[^i+^^l[a6*|(/.i)(02|+«*fo|02)(</>i|]. 

We see that the off-diagonal elements decay at a fixed rate 
while the diagonal elements remain constant. It satisfies 



d , . 



1, 



-[A1 + A2] p(t)-^Q»pWQ. 



where the projection operator Q.; 



that this is a Lindblad equation, note that Q 



To see 

= Q^ 



and J2i=i QlQi — 1- There are two Lindblad operators: 
identify V = ^(ATTA^Q. in F.q.Q. 

B. Unitary Jump 

Suppose in time dt, a state vector has probabil- 

ity Xdt of changing to exp —iG\ip, t) (probability 1 — Xdt 
of being unchanged), where G is a hermitian operator 
and exp— iG ^ 1. The density matrix at time t -I- is 
therefore 

p{t + dt) = (1 - Xdt)p{t) + Xdte-'^ p{t)e''^ , 
so its evolution equation is 
d 



dt 



p{t) = -X[p{t) ~ e-^^p{ty 



This is of the Lindblad form, with one Lindblad operator 
L = ■\/2Aexp— iG. 

In the basis where G is diagonal with elements (ffi, 52), 
we get dpii/dt = and 



dt 



Puit) = -Api2(t) 



1 



So, again, its diagonal elements remain constant. Its off- 
diagonal elements decay at the fixed rate A[l — cos(52 — 
51)] and their phases change. 



C. Random Unitary Transformation 

Suppose in time dt a state vector \tp,t) undergoes a 
unitary transformation to 



1 



\^,t) = [l-tG9--iG9y + ...\4,,t) 



with probability 

P{e)dd 



de 



y/AirXdt 
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The density matrix a,t t + dt (neglecting terms of order 
higher than dt) is given by 



p{t + dt) 



P{0)d9e-'^'^p{t)e 



iG9 



Pit) 



p{9)de 

Xdt 



p{t) - [^G2p(<) + p[t)\G^ - Gp{t)G 



G'p{t)+p{t)G^ -2Gp{t)G 



giving the Lindblad equation 



dt 



Pit) 



G^pit)+p{t)G^-2Gp{t)G 



[G[G,p(t)] 



with one Lindblad operator L = vAG. 

In the basis where G is diagonal with elements (ffi, 52), 

^Pvit) = -^iSi - 9])^Pijit)- 

So, again, its diagonal elements remain constant but its 
off-diagonal elements decay at a rate determined by the 
difference in eigenvalues. 



D. State Exchange 

Suppose, in time dt, with probability Xdt, a state vec- 
tor \^j},t) exchanges its basis states ,\<f)2) , becoming 
a\il),t), where {<pi\a\4>j) is the Pauli matrix with diagonal 
elements and off-diagonal elements 1. It is easy to see 
that the density matrix evolution equation is 
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d 



pit) = -X p{t) - ap{t)a 



and is of the Lindblad form, with one Lindblad operator 
L = V2Aa. 

The density matrix elements therefore satisfy 



d , d 
d , d 

:Pl2[t) ^ -T7P2l(<) 



-X[pn{t) - P22{t)], 

-A[pi2(t)-P2l(<)]. 



dt ^ ' dt' 
The diagonal density matrix elements change in this ex- 
ample, decaying to 1/2. The off-diagonal matrix elements 
keep their real parts while the imaginary parts decay to 
0. 



E. State Transitions 

Here we consider arbitrary N . Suppose in time dt, a 
state vector makes a transition to state 

(nlV',^) 



TO) 



IHV,i)| 



with probability 



p„Adt|(?/;,i|n)|^ 



(Pm > 0, X]m=iP™ ~ ^^'i probability of all such 
transitions is X^mn Pmn = ^dt, so the state vector is un- 
changed with probability 1 — Xdt. Define Qmn = |TO)(n|. 

Note that I]m,n=l P™QL«Qmn = 1- 

The density matrix at time t + dt is 

N 

p{t + dt) = (1 - Xdt)p{t) + Xdt ^ PmQmnP(i)Qnm, 

m,n— 1 



so its evolution equation is 
d 



dt 



Pit) 



-A 



Pit) 



N 

E 

m.n—l 



This is of the Lindblad form, with (one more than 
the necessary maximum!) Lindblad operators L'"" = 

ran • 

The matrix elements of the density matrix obey 



d 

J^Prsit) 



-X 



-X 



Prs (t) - PrSrs ^ Prm (0 
n 

Prsit) -prSrs^ ■ 



The off-diagonal elements decay at a uniform rate. The 
diagonal elements do not remain constant. They decay 
to predetermined values Pr'. 



Prrit) = Prr{0)e 



-At 



Pr 



1 



-At 



This might be useful in modeling the approach to thermal 
equilibrium, where the states |to) are energy eigenstates 
and Pr is the Boltzmann probability exp—Er/kT. 



III. APPLICATION OF CONSTRAINTS 

We now turn to deriving the Lindblad equation as the 
most general equation satisfying the constraints. 

While the hilbert space discussed here shall be assumed 
of dimension N, N may be allowed to go to to infinity 
and, also, the argument may readily be extended to a 
continuum basis. 

To make the argument easier to follow, examples of 
how its steps apply to a two-dimensional hilbert space 
shall occasionally be inserted. 

The Markov constraint is that the density matrix 
p{t') = p' at a later time t' , depends only upon the den- 
sity matrix p(t) = p at an earlier time t, not upon the 
density matrix over a range of earlier times. 

The linearity constraint, combined with the Markov 
constraint, is that the matrix elements of p' can be writ- 
ten as the sum of constants multiplying the matrix el- 
ements of p rather than, say, powers of the matrix ele- 
ments of p or any other kind of function of these matrix 
elements: 



AT 

^ ^ ^ir^jsPrs- 
r,s— 1 



(4) 



Here p^^ = {(j)i\p'\(j)j) with some convenient orthonor- 
mal basis and, similarly, p is expressed in the same basis. 
The constants Airjs can be functions of t' , t. There are 
constants, and each can be complex, so there are 2iV^ 
real constants involved in Eq.0 



The hermiticity constraint = p^* = p'^ , applied to 



Eq.Q, results in 



N 
r,s— 1 



[A* 



(5) 



1. Two-dimensional space: hermiticity 

Suppose we have the equation 
-Biipii -I- -B12P12 + -B21P21 + -B22P22 = TrBp = 0, 

(the Bij are constants) which holds for all possible den- 
sity matrices. Then, one can see B = as follows. 

First choose the density matrix pn = 1, with all other 
elements vanishing: thus, Bn = 0. Similarly, one shows 
B22 = 0. Next, employ the density matrix pij = 1/2, 
which results in B12 + B21 = Finally, use the density 
matrix pn = P22 = 1/2, pi2 = — P21 = i/^i which results 
in B12 — B21 — and so B12 = B21 = 0. Therefore, 
B = 

The four density matrices used here, 

2(1 + -')^ li^--'h 2(1 + -'X 2 + 

(written in terms of the Pauli matrices) we shall call the 
density matrix basis. Any 2x2 matrix can be written 
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as a linear sum with constant (complex) coefficients of 
these four matrices. More than that, they form a matrix 
basis for hermitian matrices, in that any hermitian ma- 
trix can be written as a linear sum with constant (real) 
coefficients of these four matrices. More than that, and 
this is the reason for their deployment here, they form 
a matrix basis for density matrices, in that any density 
matrix can be written as a linear sum with constant pos- 
itive real coefficients of these four matrices such that the 
sum of the coefficients add up to 1. 

This basis is to be distinguished from another basis, the 
Pauli matrices plus the identity matrix, which we shall 
call the Pauli+1 basis. This is also a matrix basis for 
hermitian matrices but it is not a density matrix basis. 



2. N-dimensional space: hermiticity 



Generalizing, if we have an equation 

N 



BsrPrs = TrBp = 



(6) 



for a matrix B, which holds for all valid p, then B = 0. 

This can be seen by using an A'^^-size density matrix 
basis, (generalizing the 2^-size density matrix basis of 
the previous section) . First choose pkk ~ 1 with all other 
elements vanishing, which implies Bkk ~ 0. Then for 
particular values of k,l, choose p^k = Pii = Pki = Pik ~ 
1/2 with all other elements vanishing, from which one 
finds Bki + Bik = 0. Finally, choose pkk = Pii = ipki = 
—ipik — 1/2, from which one finds Bki — Bik — 0, so 
Bki = Bik — 0. Letting k, I range over all possible pairs 
of indices results in B = 0. 



Two-dimensional space: evolution equation and trace 
constraint 



It therefore follows from Eq.Q that 



A* — 4 - ■ 

where each index can take on the values 1 or 2. 

A matrix B for which B^ „ — Bn^m is a hermitian 
matrix. Therefore, A is a hermitian matrix, where we 
regard the number pairs 11, 12, 21, 22 as four different 
indices. That is, A is a 4 x 4 dimensional matrix. The 
most general 4x4 hermitian matrix is characterized by 
16 real numbers (the four real diagonal matrix elements 
and the six complex matrix elements above the diagonal). 
Since there are 32 real numbers which characterized the 
most general superoperator in a two dimensional space, 
the condition of hermiticity of the density matrix has cut 
that number in half. 

A hermitian matrix can be written in terms of its or- 
thonormal eigenvectors and eigenvalues, and that decom- 
position shall prove very useful here. There are four real 



eigenvalues. A", where a = 1,2,3,4. Corresponding to 
each eigenvalue is an eigenvector E" in the four dimen- 
sional complex vector space. 

The four complex components of E^^. make 32 real 
numbers, but they are constrained. Each eigenvector 
is normahzed to 1: J^f r=i E^^, = 1 + iO provides 8 
constraints, lowering the number of free components to 
24. The orthogonality of to the other three vectors 
provides 6 constraints, the orthogonality of E^ to the 
remaining two vectors provides 4 constraints, and the or- 
thogonality of E^ to E"* provides two constraints. Thus, 
there are 12 constraints on the 24 free components, so the 
eigenvectors contain 12 free components. These, together 
with the four eigenvalues, comprise the 16 real numbers 
characterizing A. 

An example of such an orthonormal basis is given by 
l/\/2 multiplying the Pauli-|-1 basis. If we write the 
four components of E" as a four dimensional vector 
with components [E",E1'2, E^^,E^^], then E^ = 2-^/'^cr^ 
has components 2^i/2|-o^ g], E^ = has com- 

ponents 2^^/^[0, — i, z, 0], E^ = 2~^/^cr^ has compo- 
nents 2" 0,0,-1], E'' = 2-1/21 has components 
2-i/2[l, 0, 0, 1], It is easy to verify that this is an or- 
thonormal set of vectors. 

Although each E" is a vector in a four dimensional 
space, with four components E'j" , E" can also be regarded 
as an operator in the two-dimensional hilbert space with 
four matrix elements E°j. This leads to a neat way of 
writing the orthogonality relations for these eigenvectors. 
Instead oi J2iLiY^^=i Efr^i* = 5°-^ , we can write 

TrE"E^^ = 5"l^, 

where E^^ is the hermitian conjugate (complex conjugate 
transpose) of E'^. It is easy to see how this works for the 
example where E" is l/-\/2x the Pauli-|-1 basis. 

The expression for the components of A written in 
terms of its eigenvectors and eigenvalues is 



a=l 

Putting this into Eq. Q results in the evolution equat 

4 2 4 



ion 



a—l r,s— 1 a — 1 

Now, lets impose the trace constraint, i.e., Y^^=iP'ii ~ 
1. In terms of components this says 

4 2 2 

i = E^"E E E^prsE*^- 

a — l i—1 r,s— 1 

Writing 1 = X^Li = J2l,s=i ^srPrs, the trace con- 
straint can be written as 



E 

r,s—l 



a—l i—1 



Prs = 



5 



or m matrix notation as 

4 



Tr 



A"E"^E" - 1 



Q = l 



where 1 is the unit matrix. This must hold for arbitrary 
p. We have seen how to handle such an expression. By 
successively putting in the four density basis matrices, 
we obtain the trace constraint 



4.. N- dimensional space: trace constraint 

The A^-dimensional case works just like the two- 
dimensional case. It follows from Eq.([5]) that A can be 
viewed as an N"^ x iV^ hermitian matrix. It has real 
eigenvalues. Its N"^ complex eigenvectors Ef^ satisfy the 
orthonormality conditions 



IV. COMPLETE POSITIVITY 

The final constraint is positivity. This says, given an 
arbitrary N-dimensional vector \v), that the expectation 
value of the density matrix p' is non-negative. This con- 
straint, applied to Eq.(l8|), is 



< {v\p'\v) = ^ A"(w|E"pE 



where we have defined E"'''|w) = \vo 
Positivity of p ensures {va\p\va) 



E 

a=l 



Va\p\Va), 
(10) 



> 0. Thus, we see 
from Eq.(|lO|), if all the A"'s are non-negative, then p' 



will be positive too. 

However A" > 0, while just shown to be sufficient for 
p' to be positive, is not necessary. In the next section, 
we shall give an example where an eigenvalue is negative, 
yet p' is positive! 

Therefore, a stronger condition than positivity is nec- 
essary to ensure that A" > 0. This condition, presented 
after the example, is complete positivity. 



N N 



(7) 



i—l r—1 



With A written in terms of its eigenvalues and eigenvec- 
tors, Eq.Q becomes the evolution equation 

N 

P'^3 = E E or 
a— 1 r,s—l 

p' = ^A^E^pE"!'. (8) 

Q=l 

A" and E" depend upon t' — t, but we shall not write 
that dependence until it is needed. 

Next, imposition of the trace constraint on Eq.Q, with 
Trp' = 1 = Trip, gives 



5. Two-dimensional space: example of a positive density 
matrix with negative eigenvalue 

This example uses the Pauli+1 eigenvectors E" ~ 
2"i/2^Q and 2~^^^1. (Note that the trace constraint 
^ is satisfied, provided J2a=i ~ 2' since the square 
of each of the Pauli-l-1 matrices is 2~^/^l.) Choose 

Ai = A2 = -A3 = A4 = 1: 



Ipl] = 



P22 P12 
P21 Pii 



p' is just p with its diagonal elements exchanged. Thus, 
because p is positive, then p' is positive. This is a partic- 
ularly simple example of a more general case discussed 
in Appendix [S] 



Tr [^A"E"^E" - 1 



Using the density matrix basis as in Eq.([6| et seq., we 
obtain the trace constraint: 



^ A"E"^E" = 1. 



(9) 



By taking the trace of Eq. (|9| and using Eq.0, we find 
the interesting relation 



^ A" = iV. 



6. N-dimensional space: definition of complete positivity 

It is not positivity but, rather, complete positivity that 
makes the non-negative eigenvalue condition necessary. 
Here is what it means. 

Add to our system a non-interacting and non-evolving 
additional system in its own iV-dimensional hilbert space. 
The enlarged hilbert space is of dimension N"^ . The sim- 
plest state vector in the enlarged space is a direct product 
\(j)i)\Xj)'- \4'i) is a vector from the original hilbert space, 
\Xj) is a vector from the added system. The general state 
vector in the joint space is the sum of such products with 
c-number coefficients. 

Form an arbitrary density matrix TZ for the enlarged 
system. Suppose it evolves according to Eq.Q, where 
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E" is replaced by E" x 1 (i.e., the evolution has no effect 
on the vectors of the added system.) Complete positivity 
says that the resulting density matrix TV must be posi- 
tive. 



7. Two-dimensional space: complete positivity 

Complete positivity says, given the evolution equation 
([8|, that {w\'RJ\w) > for an arbitrary dimensional 
vector |w) and for any initial density matrix TZ in the 
enlarged hilbert space. We wish to prove that complete 
positivity implies the eigenvalues are non-negative. What 
we shall do is judiciously choose a single vector and 
four pure density matrices TZ so that the expressions 



Qfl 



•\Tl'\w) = ^ A"('u;|E"7^E 



are ^ A*^, with a positive constant of proportionality. 
Therefore, for complete positivity to hold, A*^ must be 
non-zero. Here are choices that will do the job. 
We shall choose the maximally entangled vector 



El 



{{w\'w) = 4, but it need not be normalized to 1). We 
construct the state vectors 

and use them to make four pure density matrices 
\^)'^ ^ii^]. (Note that Tr^)l^'^{^\ = 1 because of the 
orthogonality relation Eq. (?])). Then, for one /3, 

4 

Putting this into Eq.([8]), the complete positivity condi- 
tion is 

4 

< {w\TZ'\w) = J2 
(0,|(x.|E"|0,)|x,)(0d(x,'|E^"l0r')IXr') 

a,i,j^i' ,j' — 1 
4 

= A"Tr(E^^E")Tr(E'3E""f) 

a=l 

4 



^A"(5"^)2 = A^ 



a=l 



(using the orthogonality relation ([7])). 
Thus, complete positivity implies A*^ > 0. 



8. N-dimensional space: complete positivity 

We follow the same procedure in the N-dimensional 
case. However, to be a bit more general, we shall use 
an arbitrary vector jw), and an arbitrary pure density 
matrix TZ: 



N 



m,n— 1 
TV 



(lla) 



7^ ^ CkiC*,,M\xi){M{xi'\ (lib) 

k,l,k',l' = l 

where Cki , Dmn are yet to be spec ified complex con- 
stants. The unit trace of p in Eq.( lib ) requires TrC^C = 



1. Then, the complete positivity condition is 
< (^^;|7^V) ^ ^ A"(^^;|E"7^E"V) 

AT^ N 

a=l 1 

= ^ A"rr[CD^E"]Tr[E"^DC^]. (12) 



Now, choose DC''' = E^, for any particular /3. This 
choice can be made in many ways. Two are C^^ — , 
D = 1 (the choice made in the two-dimensional exam- 
ple just discussed) or = E'^, C = N'^^'^l (note, both 
choices respect Tr C^C — 1 ) . With this choice in Eq. ( 1 2 1 , 
and with use of the orthonormality conditions Eq. Q , we 
obtain as the consequence of complete positivity: 



< ^ A"((5"^)2 = A*^ for 1 < /3 < 7V2_ (^3) 



V. KRAUS REPRESENTATION 

We have now applied all the constraints needed to ob- 
tain a valid density matrix p' at a later time t' from an 
earlier density matrix p at time t. This relation is Eq.(|8|, 
supplemented by the orthonormality conditions Q, the 
trace constraint ^ and the condition of non-negative 
eigenvalues (13). 

It is customary to define M" = %/A"E", so that Eqs.^ 
|9l) can be written in terms of M" alone: 



9' ^Y M"pM"1' 



^M"'fM" = 1. 



(14a) 



(14b) 



(However, the orthonormality conditions, written in 
terms of M", now depend upon A"). Eq.(14a) is called 
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the Kraus representation and {M"} are called Kraus 
operators[5]. 

We have proved the necessity of the Kraus represen- 
tation, but it i s also su fficient. That is, for any {M"} 
satisfying Eqs.( 14a|l4b ), even for more than iV^ opera- 
tors, also with no orthonormality conditions imposed, all 
the constraints on p' are satisfied. It is easy to see that 
hermiticity, trace 1 and positivity are satisfied. Complete 
positivity requires a bit more work, and that is given in 
Appendix [B] 

This general statement of the Kraus representation 
might seem to imply a larger class than we have derived 
as necessary, but that is not so. Since the Kraus repre- 
sentation is hermitian, trace 1 and completely positive, 
it may be written in the form Eq.([8|), as we have shown. 



For /? 7^ iV^, the eigenvalues vanish. Note that the 
condition TrE^l = says that these eigenvectors are 
orthogonal to ^ 1. 

And, indeed, in this case, Eq.(|8| becomes the identity 



p{t)^lp{t)l = N^lpit)^l. (18) 



B. Eigenvectors and eigenvalues when t' — t + dt 

When t' = t + dt, the eigenvalues and eigenvectors 
change infinitesimally. Accordingly we write 



VI. LINDBLAD EQUATION 

Now that we have satisfied all the constraints on the 
density matrix p' = p(t'), we can let t' ~ t + dt, and 
obtain the differential equation satisfied by p{t). For the 
rest of this paper we shall only treat the A'^-dimensional 
case since the argument is precisely identical for the two- 
dimensional case, except that N = 2. 



A. Eigenvectors and eigenvalues when t' = t 

First, lets see what we can say about the eigenvectors 
and eigenvalues when t' = t. Then, Eq.([8| says 



N 



Pn 



^Y.^"I1 ETrK^Prs or 



a— 1 r,s — l 



N 

E 



a=l 



(15) 



As we have done before, successive replacement of p by 
the N"^ members of the density matrix basis results in 



(16) 



Multiply Eq.(16) by E^^ and sum over 7, s Use of the 
orthonormality relation ([7| gives 



(17) 



If TrE/^ ^ and A'' ^ 0, Eq.(17) says that all the 



eigenvectors arc ^ 1. But only one of a set of orthog- 
onal eigenvectors can be proportional to the identity. 
Therefore, for the rest of the eigenvectors, A^ = and 
Tr^P = 0. 



Can one eigenvector E^ = N"^l'^l. From Eq.(ll7j), 
we find the associated eigenvalue A^ = N . 



A^ [dt] = N[l - dt], X°'{dt) = c°'dt {a ^ N^), 

E^\dt) = ^[1 + Bdt], E"(di) ^K"{ay^N^), 
V N 



(19) 



where the c" are constants. We do not include a term 
^ dt in the expression for E"{dt) since, because X°'{dt) ~ 
dt, i t w ould contribute a negligible term ~ (dt)"^ to 
Eqs.([8l9l). 

Because the eigenvalues must be positive, and because 
the eigenvalues sum to N (equation following Eq.([9])), we 
see that c" > (all a). B and K" are restricted by the 
orthonormality conditions, which we shall look at later. 



C. The evolution equation 



Putting Eqs.(19) into the evolution equation iSl gives 



p{t + dt) = [1 - cJ^ dt][l + Bdt]p{t)[l + BMt] 
+dt c"K"p(i)K"^ or in the limit dt 0, 



a=l 



dt 



Pit) 



' p{t) + Bp{t)+ p{t)B^ + c"K"p(t)K"t. 

(20) 

Putting Eqs.(19l into the trace constraint ([9| gives 

A^[l - c'^^dt] [1 -f BMt] [1 -f Bdt] 

+dt c"K"^K" = 1, or, 



1 = B + Bt + Y c"K"tK". 



(21) 
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Using (21) to replace c in (20 1 (specifically, c p 



(l/2)[c^ Ip + pc^ 1]) results 



d 
dt 



1 ^ 

- J2 c"[K"K"VW + p(i)K"K"^ - K"p(t)K"^]. 

(22) 



a=l 



D. The Lindblad equation 

If we define -iU = (1/2) (Rt - B) and L" = Vc"K", 
the evolution equation ( 22 1 becomes the Lindblad equa- 
tion ([3]): 



dt 



Pit) 



H,p(t) 



N^-l 



-- [L"L"V(i) + p(OL"L"^ - 2L"p(t)L"^ 



(23) 



E. Orthonormality conditions 



It is a consequence of this derivation that the Lind- 
blad operators L" in Eq.(23) are not arbitrary opera- 



tors, because they are restricted by the orthonormality 
conditions ([7|. Putting Eqs.(19) into Eq.([7| constrains 
B, K". 

For what follows, we recall from the discussion in sec- 
tion II that {K"} can be regarded in two ways. In one 
way, they are regarded as N'^ — 1 operators acting on 
vectors in an N dimensional space, with matrix elements 
^i^j ihj — 1---/V)- In the other way, they are regarded 
as — 1 vectors in an iV^ dimensional space, each with 
components {Ki^, K12, ■■■K'^j^). In particular, the trace 
of two operators is the same as the scalar product of two 
vectors, as in Eq.([7]) 

The orthonormality relation ([?]), applied successively 
to (a = (3 = N^)M ^N^p^ N'), {a ^ p ^ N^), 



with use of Eqs.(19l, are 



Tr[B + B^]=0, (24a) 
TrK" 0,, (a = l,...iV2 - 1) (24b) 
TrK^K'^l' = (5"'3 (a, P ^ I, ...N^ - 1). (24c) 



Eq.(24a) says that the hermitian part of B vanishes. 



This provides no restriction at all on H, which is the 
anti-hermitian part of B. 

Eq.( 24c[) says that the vectors K" are orthonormal. 

Eq.(24b|) says that TrK"! = 0, which implies that 
N~^^^l completes the orthonormal set. 



F. The general Lindblad form 



We shall now show that the Lindblad equation (23) 



with arbitrary Lindblad operators (no constraints what- 
soever) can be transformed to new, constrained, Lind- 
blad operators of Eq. ( |22| ) by adding a constant (achiev- 
ing the vanishing trace constraint (24b)) followed by a 
unitary transformation(achieving the orthogonality con- 
straint ( |24c[ )). 

First, we see that we can transform the arbitrary Lind- 
blad operators L" to Lindblad operators L'" which are 
traceless. Define L" = L " + where the fc" are 

iV^ — 1 constants, and substitute that into the Lindblad 
equation, obtaining : 



dt 



Pit) 



U + i{k"L'"'< - fc"*L'"),p(t) 



1 ^ 

-- ^[L'"L'"V(i) + p(i)L'"L'"^ - 2L'"p(t)L'"^]. 



With a redefinition of H, this is again the Lindblad 
equation, expressed in terms of L'". By choosing fc" = 
iV~^TrL", the new Lindblad operators satisfy TrL'" = 
0. 

Now write the — 1 Lindblad operators L" (hereafter 
assumed traceless) in terms of — 1 new operators 
(a, /3 = 1, ...N'^ — 1) using the linear transformation 



E 

/3=1 



(25) 



We ask what the matrix [/"'^ must be in order that 
the Lindblad form be unchanged (L" replacing L" in 
Eq.(23)). With arbitrary operators Ai, A2, A3, 



J2 AiL"A2L"U3 = U"^^U*'''^'Ail^A2t^'^A: 

a,l3,P' = l 



3- 



(The three terms in the Lindblad equation have two 
of the A = 1 while the third A = p). This equals 

^^^7^ AiL'^A2L'^^A3, leaving the Lindblad equation 
unchanged in form, if and only if the matrix C/"'^ is uni- 

Inverting Eq.(25l, we see that is traceless. 



Now, consider the matrix TrL"L" It is hermitian, so 
its eigenvalues are real and its eigenvectors are orthog- 
onal. It can be brought to diagonal form by properly 
choosing our unitary transformation, so we obtain 



Trt^L'^"' = y [/t'3."[/"'./'rrL"L"'^ = (26) 



This is almost the orthogonality constraint (24c). The 
eigenvalues are non-negative, since it follows from 
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Eq.(26l that 



N 



Thus, according to Eq.(26), the L'^ are orthogonal vec- 
tors, with squared norm c*^ . We can add one more vector 
^ 1 to complete the set, orthogonal to the rest since 
Trt^l = TrLf^ = 0. 

We can define new operators K*^ which are orthonor- 
mal and traceless, by = V&^K-^. In terms of these 
operators, the Lindblad equation (23) written in terms 
of Li^ becomes 



dt 



Pit) 



H,p(t) 



(27) 



a = l 



This is precisely Eq.(22| with K", c", replacing K°, c° . 
Moreover, the orthonormality constraints Eqs.( 24b|24c ) 
on K° are satisfied by K". 



G. Concluding remarks 

We have shown that the constraints on the density 
matrix p' mandate its evolution equation (22), where 
c" > and B,K" satisfy the constraints (24). We then 



showed that this is completely equivalent to the Lindblad 
Eq. ( 23 1 , with no constraints at all on the N'^ — 1 Lindblad 



operators L". 

However, there is no need to restrict the Lindblad 
equation to no more than N'^ — 1 operators. We conclude 
our presentation by showing that the Lindblad equation 
with any number of operators has all the required prop- 
erties. (Of course, from what we have shown, such an 
equation may be reduced to one with no more than iV^ — 1 
operators.) 

Looking at the Lindblad equation in that case, 
j^p{t) = ~t\ii,p{t) 

-\ 5^[L"L"tp(t) + p(i)L"L"t _ 2L>(t)L"t](28) 

OL 

it is easy to see that hermiticity, and trace 1 (in the 
form dTrp{t)/dt = 0, with Trp{to) — 1) are satisfied. 
Complete positivity requires a bit more work. 

We have shown in Appendi x [B] t hat complete positiv- 
ity holds for the Kraus form (14a) subject to the trace 



constraint (14b), with an arbitrary number of operators. 
So, if Eq. ( 28 ) can be written in the Kraus form with the 



Accordingly, we write Eq.(|28[) as 

1 



p{t + dt) = [l - dt{m + 2 H L"^L")J p{t) 

a 

■ l-dt(-iH-f i^L"^L")] +dt^L"p(i)L"^ 

a a 

(29) 

Identifying the Kraus operators as 



M° = [1 - dt{m +1^ L"^L") , M"^" = VdtL", 



gives the Kraus form Eq.( |14a[ ), 

p{t + dt) = ^M"p(t)M"^ 



(30) 
(31) 



Now, take the trace of Eq.(31). Since Eq.(28) implies 



Trp{t + dt) = 1 = Trlp{t), the resuh is 
1 = rr^M"^M"/9(t) or Tr [^M"^M" - 1 p{t). 

a a 

(32) 

As we have done before, by successively replacing p(t) by 
the members of the density matrix basis, we obtain the 
Kraus trace constraint ( |14b I: 



(33) 



Therefore, the Lindblad form with an arbitrary num- 
ber of Lindblad operators, is completely positive. 



Appendix A: A Class of Positive Density Matrices 
With Negative Superoperator Eigenvalues 

Consider the Pauli+1 matrices (multiplied by l/\/2) 
as the eigenvectors E^. As noted in Section IV, they 
satisfy the trace constraint ^ if J2a=i ^ 2. Eq.([8|) 
becomes 



p' = Y A^E^pE^I' 



a=l 



-[AVVct^ + AV^pcr^ ^ X3fj3pfj3 + A^lpl] 

(Ai + A2)p22 + (A3 + A4)pn, (Ai - X^)p2i ~ (A^ - X^)p,2 
(Ai - A2)pi2 - (A3 - A4)p2i, (Ai + A2)pn + (A^ + X'^)p,2 

(Al) 



trace constraint, we have shown it is completely positive. 



Now, we can successively replace p by the four density 
basis matrices, and demand that the A" be chosen so p' 
is positive for all. Since the sum of these four p's with 
positive coefficients (adding up to 1) is the most general 
two-dimensional density matrix, then the most general p' 
will be positive. 
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The first density basis matrix has pn = 1 and the rest 
of the matrix elements vanishing. Then, 



(A3 + A4) 
(A1 + A2) 



(A2) 



Since the eigenvalues of the density matrix must lie be- 
tween 1 and 0, we obtain the two conditions: 



2 > A^ + A^ > 
2 > A^ + A"^ > 0. 



(A3a) 
(A3b) 



The second density basis matrix, with p22 = 1 and the 
rest of the matrix elements vanishing, gives the same re- 
sults. 

The third density basis matrix is p — (1/2)[1 + ai]. 
Using J2a=i A" = 2 to simplify the result, we obtain: 



(Ai + A^-l) 



(Ai + A* - 1) 
1 



(A4) 



The eigenvalues of p' here are {X^+\'^)/2, l-(Ai + A'^)/2, 
so the condition that they lie between and 1 is 



2 > A^ + A"^ > 0. 



(A5) 



The fourth density basis matrix is p = (1/2)[1 + a-2\. 

4 



Using 'Ylia=\ A" = 2 to simplify the result, we obtain: 

1 



P 



1 z(A2 ^ A^ 

^z(A2+A4-1) 1 



(A6) 



The eigenvalues of p' here are (A2 + A4)/2, 1-(A2 + A'^)/2, 
so the condition that they lie between and 1 is 



2 > A^ + A"^ > 0. 



(A7) 



So, we have obtained the result that p' will be positive 
if Eqs.( |A3a|[A3b|A5|A7[ ) and 



a=l 



(A8) 



are satisfied. The sum of Eqs.(A3b, A5 A7))minus (A8) 
tells us that 2 > A"^ 



> -1. 



Eq.(A8) and the constraint boundaries are three- 



dimensional hyperplanes in the four dimensional A-space. 
Their intersections delineate the allowed areas for the 
eigenvalues. 

We shall be content here to set A^ = A^ (in which case 
Eqs.(A3a) simphfies to 1 > A^ > 0). Then, Eq.(|A8| de- 
scribes a plane in {A^, A'^, A*} space, and its intersection 
with the constraint boundary planes can be drawn. This 
is shown in Fig. 1. There are two regions where one of 
the eigenvalues is negative and the other two are positive: 
the points in the heavily outlined upper left triangle have 
A^ < 0, and the points in the heavily outlined lower right 
triangle have A^ < 0. 



A' =2 




A* =-1 



FIG. 1: Allowed regions of eigenvalues for a positive den- 
sity matrix p' . The two dark-outlined triangular regions are 
where an eigenvalue is negative. They abut an isosceles tri- 
angle, the restricted region of complete positivity, where all 
the eigenvalues are positive. 



Appendix B: Complete Positivity of the Kraus 
representation 



The Kraus form Eq.(14a) and the Kraus constraint 



Eq. ( 14b I , generalized to any number of arbitrary opera- 
tors M", are respectively 



f) ^ M"pM"^ 

a 



(Bla) 
(Bib) 



We want to show complete positivity. Call any one of the 
M" = M. If we can show complete positivity for MpM^^ 
for an arbitrary M, then Eq.(Bla), which involves a sum 



of such terms, will be completely positive. And, it is only 
necessary to prove complete positivity for 7^ = p x p^, 
where p^ is any possible basis density matrix (described 
in the paragraph following Eq.([6])) in the added hilbert 
space, since the most general density matrix in the direct 
product hilbert space is the linear sum of such terms. 
We now calculate 



{v\VJ\v) = (t;|MpM^ X p*^ 



(B2) 



for arbitrary \v) = Y.m,n=iDrnn\4>m)\Xn) ■ There is no 
loss of generality if we pick the basis vectors \xn) in the 
added hilbert space any way we like. We shall pick them 
to be the eigenstates of p'^. Now, we note that each p'^ 
has one eigenvalue 1 and the remaining A'^ — 1 eigenvalues 
are 0. Call the eigenvector \xi) which corresponds to the 
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eigenvalue 1. Then, 

m,m' = l 

N N 

= [ ^ D*^,,{cl)m'\M]p[ ^ DmiM^ld^m}] > 0, (B3) 



N 



m'=l 



m'=l 



Therefore, the Kraus form with arbitrary M"'s is com- 
pletely positive. 
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